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ABSTRACT 

The X-Ray Telescope (XRT) on the Japanese/USA/UK Hinode (Solar-B) 
spacecraft has detected emission from a quiescent active region core that is con- 
sistent with nanoflare heating. The fluxes from 10 broadband X-ray filters and fil- 
ter combinations were used to constructed Differential Emission Measure (DEM) 
curves. In addition to the expected active region peak at Log T = 6.3-6.5, we 
find a high-temperature component with significant emission measure at Log T 
> 7.0. This emission measure is weak compared to the main peak - the DEM is 
down by almost three orders of magnitude - which accounts of the fact that it 
has not been observed with earlier instruments. It is also consistent with spectra 
of quiescent active regions: no Fe XIX lines are observed in a CHIANTI syn- 
thetic spectrum generated using the XRT DEM distribution. The DEM result is 
successfully reproduced with a simple two-component nanoflare model. 

Subject headings: Sun: corona, Sun: X-ray radiation, Sun: fundamental param- 
eters 



1. Introduction 

If nanoflares heat the corona (Parker 1983, 1988; Klimchuk 2006), then we would expect 
an ever-present hot component with a temperature of 10 MK or higher, even in quiescent 
active regions. A glimpse of this hot material may have been seen by, for example, the Hard 
X-Ray Imaging Spectrometer on Solar Maximum Mission (Martens, van den Oord & Hoyng 
1985) or the Mg XII imager onboard CORONAS-F (Grechnev et al. 2006). For the most 
part, however, it has only been possible to put an observational upper limit on the amount 
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of this hot material because high temperature spectral lines like Ca XIX (Lemen, Bentley & 
Sylwester 1986; Fludra & Schmelz 1999) and Fe XIX (Schmelz 1993; Wang, Innes & Solanki 
2006; Teriaca et al. 2006) were not seen except during solar flares. Even observations of the 
pervasive S XV line by the Yohkoh Bragg Crystal Spectrometer (Culhane et al. 1991) could 
be explained with cooler (3-4 MK) plasma (Watanabe et al. 1995; Sterling 1997). 

Broadband imagers like the Yohkoh Soft X-Ray Telescope (Tsuneta et al. 1991) were 
designed to investigate hot, flaring plasma, but results tended to rely on filter ratios and an 
isothermal approximation to estimate plasma parameters. In this letter, we use data from 
the X-Ray Telescope (XRT; Golub et al. 2007) on Hinode (Kosugi et al. 2007). XRT is a 
broadband instrument similar to the Soft X-Ray Telescope, but with better spatial resolution, 
sensitivity, and temperature coverage. A detailed 3 MK thermal map of an active region 
observed by XRT has already been published by Reale et al. (2007) using an innovative 
combined filter ratio method. The properties of the instrument allow us to go one step 
further, however, and perform a true multithermal analysis. XRT should be able to observe 
the high temperature plasma predicted by nanoflare models, or at least determine a strict 
upper limit for the amount of emitting material. 



2. Analysis 

XRT observed AR 10955 (S09W30) on 2007 May 13 at 1800 UT. The data set includes 
full Sun, full resolution images in all filters as well as two filter combinations. The data were 
processed with standard XRT software available in SolarSoft (xrt_prep.pro). This includes 
subtraction of a model dark/read- noise frame, correction for vignetting, removal of high 
frequency pattern noise, and normalization by exposure time. Spacecraft jitter was removed 
(xrt_jitter.pro), and long and short exposures (see Table 1) were co-aligned and combined to 
increase image dynamic range. The data from the thickest channel (Be_thick) were subject 
to additional Fourier filtering to remove some low-level, residual, longer-wavelength noise 
patterns. We have used updated filter calibrations (Narukage et al. 2009) and 1640 A of 
diethylhexyl phthalate, a time- dependent contamination layer on the CCD, which can now 
be modeled using the make_xrt_wave_resp.pro and make_xrt_temp_resp.pro routines available 
in SolarSoft. 

Figure 1 shows the XRT TLpoly image of AR 10955, which produced a small GOES A7 
flare at 18:12 UT. Our observations indicate that the effects of this flare were limited to the 
left side of the active region, so we focused our attention on the outlined area. Even here, well 
away from the flaring material, there are significant counts in the Be_thick channel. Since 
we are specifically searching for hot plasma, we restricted our analysis to those pixels within 
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the outlined area having significant counts in Be_thick. The noise in the Be_thick image 
(on the disk but outside the active region and bright point areas) is: — 0.00039±0. 00245 
Data Numbers s _1 and we chose a threshold of > 0.03 Data Numbers s _1 . Given the known 
response of XRT to optically thin thermal plasma, we can generate the possible Differential 
Emission Measure (DEM) curves that can reproduce the observed fluxes in 10 different XRT 
filters and filter combinations. We chose uncertainties of 3%, which result from comparing 
the Monte-Carlo realizations of the full area outlined in Figure 1 with the individual pixel- 
by-pixel results: the 1 a spread in the two DEM distributions is similar. (This result will be 
examined in greater detail in a future paper.) 

The orange (heavy dashed) curve in Figure 2a shows the optimal median of the Monte- 
Carlo iterations generated by the Markov-Chain Monte Carlo (MCMC) based reconstruction 
algorithm (Kashyap & Drake 1998). The optimal median is the iteration best matching the 
median DEM solution in each temperature bin, after discarding non-convergent solutions. 
The routine fits a locally smoothed DEM curve to the data by comparing the predicted to 
the observed fluxes and modifying the solution randomly to obtain new realizations. The 
smoothing scale varies with temperature to account for changes in the information content 
available as codified in the filter response curves. The new realizations are kept or discarded 
according to the Metropolis criterion based on changes in the x 2 values at each step, resulting 
in a Markov-Chain that explores the parameter space efficiently. This results in a sample of 
solutions that is representative of the actual probability distribution of the DEM (see e.g., 
Smith & Roberts 1993; Casella 1996). Many of the details and advantages of this method, 
including DEM dynamic range, uncertainties, and smoothing, were discussed by Schmelz, 
Kashyap & Weber (2007) and references therein. The predicted-to-observed flux ratio for 
each filter is listed in Table 1. 

The yellow (heavy solid) curve in Figure 2a is generated with the same XRT data, 
but using xrt_dem_iterative2.pro. Testing and validation of the method with synthetic data 
leading up to the launch of Hinode were presented by Golub et al. (2004) and Weber et al. 
(2004). The routine employs a forward- fitting approach where a DEM is guessed and folded 
through each response to generate predicted fluxes. This process is iterated to reduce the x 2 
between the predicted and observed fluxes. The DEM function is interpolated using several 
spline points, which are directly manipulated by mpfit.pro, a well-known and much-tested 
IDL routine that performs a Levenberg-Marquardt least-squares minimization (Levenberg 
1944; Marquardt 1963). There are Nj - 1 splines, representing the degrees of freedom for 
Nj observations. This routine uses Monte-Carlo iterations to estimate errors on the DEM 
solution (grey dashes in Figure 2a). For each iteration, the observed flux in each filter was 
varied randomly and the program was run again with the new values. The distribution of 
these variations was Gaussian with a centroid equal to the observed flux and a width equal to 
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the uncertainty (3%). The predicted-to-observed flux ratios are also listed in Table 1. We are 
pleased to find that the two methods are in excellent agreement: the xrt_dem_iterative2.pro 
results are within la of the MCMC results for five of the XRT filters and within 2a for the 
remaining filters. 

To evaluate the quality of DEM reconstruction that can be achieved with XRT, we have 
compared the results of xrt_dem_iterative2.pro with a known input model. In each case, a 
nominal observation is calculated for a given DEM model, and then the procedure samples 
the observation 100 times by including random photon noise at the 3% level. The distribution 
of calculated DEM curves (relative to the known DEM model) indicates the accuracy and 
robustness of the analysis method. Figure 3 (top) shows the results for the isothermal cases: 
Log T = 6.1 and 6.4. Each plot shows the input delta function (solid line), the 100 Monte- 
Carlo realizations (grayscale), and the median of all realizations for each temperature bin 
(diamonds). The results use all nine filters and demonstrate that the routine can determine 
the correct result by placing most of the power in the temperature bin corresponding to the 
delta function model, with a smaller amount placed in the adjacent bins. Recall that the 
fitted DEM is a spline curve through several knots; although this is not the best function 
for fitting a delta function, our results are still quite good. For a multi-thermal example, 
we analyzed the active region model from CHIANTI (v. 5. 2.1; Dere et al. 1997; Landi et 
al. 2006). The results suggest how well the input DEM can be reconstructed as a function 
of the number of observing channels used. Figure 3 (bottom) shows the model AR DEM 
(solid line with two humps). The first panel uses only the four thickset XRT filters. It can 
determine the presence of the hotter peak, as indicated by the convergence of the median fit 
to the model DEM curve, but fails to detect the cooler material. The last panel uses eight 
filters. The cooler model component is detected, but not accurately represented, which is 
consistent with the results from Figure 2a where the Monte-Carlo iterations begin to diverge 
at low temperatures. 

The DEM realizations seen in Figure 2a are constrained by 10 XRT filters and filter 
combinations available for the active region core. The resulting curves show the expected 
component at Log T = 6.3-6.5, but there is also a high-temperature component associated 
with this emission that is down by almost three orders of magnitude from the main peak. 
This component is exactly the type of feature that we would expect if the active region were 
heated by nanoflares. The lower emission measure explains why it has not been detected by 
other instruments - the dynamic range of a typical DEM curve is approximately two orders 
of magnitude. Any features beyond this range would, in general, not be trusted. With 
XRT and its high sensitivity to hot plasma, however, this restriction does not apply. The 
peak filter responses span almost four orders of magnitude, so the DEM dynamic range is 
substantially greater than previous instruments. One important point, however, is that the 
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observations do not require that the DEM curve be bimodal. There is a possibility (see grey 
Monte-Carlo iterations) that the DEM decreases monotonically beyond Log T = 6.4. This 
has implications for the possible distribution of nanoflare energies (see below). We also note 
that all DEM results depend heavily on instrument calibration. A lot of effort has gone 
into determining the XRT filter thicknesses (Narukage et al. 2009), but problems may still 
remain. One of the advantages of DEM over traditional filter ratio methods is that these 
results are not as sensitive to an uncertainty in one particular filter. For example, both 
DEM methods used here indicate that the predicted flux for the Thick_Al filter is too high, 
a problem that could be related to the adopted thickness of the filter. Small changes to the 
filter thickness change the amount of predicted flux (and therefore the goodness of fit) but 
not the DEM shape. 

One obvious concern is that the high temperature plasma revealed in these DEM curves 
might result from the GOES A7 flare from the left side of AR 10955. The thickest filters, 
however, were observed well before the flare (see Table 1), thus reducing the likelihood of 
significant contamination. To address any residual concern, we have generated synthetic 
spectra using the ch_ss.pro routine available in CHIANTI. We have chosen the wavelength 
range of 13-14 A b ecause it contains several high temperature iron lines that can be used as 
a discriminator between flaring and quiescent plasma. The yellow DEM curve from Figure 
2a (the orange curve gets similar results) was used to produce the dotted spectrum seen in 
Figure 4. The Ne IX triplet lines (13.45, 13.55, 13.70 A) with a peak formation temperature 
of Log T = 6.2 as well as an Fe XVII line (13.82 A; Log T = 6.6) are clearly visible. These 
lines are observed routinely from quiescent active regions (Schmelz et al. 1996), but there is 
no hint of the hotter Fe XIX line (13.52 A; Log T = 6.9). The solid spectrum was generated 
with a DEM from the decay phase of a GOES B8 flare (Reeves et al. 2009). The strongest 
Fe XIX line is now obvious, as are several smaller Fe XIX and Fe XX flare lines which are all 
missing from the active region core spectrum. These results indicate that not only is the hot 
DEM component observed by XRT consistent with earlier observations from spectrometers, 
but more importantly, it is unlikely to have resulted from the flare. (Note: direct comparison 
of the AR core and A7 flare spectra would have been optimal, but lack of the thick filter 
observations during the flare prevented accurate flare DEM analysis; see Table 1 for the 
exact times.) 

We have modeled the observed DEM curve with nanoflares using the EBTEL hydrody- 
namic simulation code (Klimchuk, Patsourakos & Cargill 2008). Our goal at this time is not 
to accurately reproduce the details of the curve, but rather to demonstrate that impulsive 
heating can explain the overall shape. We therefore make a number of simplifying assump- 
tions. First, we assume the observed volume outlined in Figure 1 is comprised of unresolved 
loop strands with a half length of 60 Mm. Second, we assume that there are just two types of 
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nanoflares. Most of the strands are heated by weak nanoflares that repeat every 4.0 x 10 3 s, 
while 0.3% of the strands are heated by much stronger nanoflares that repeat every 3.0 x 10 5 
s. The model DEM curve is obtained by time averaging the weak and strong nanoflare sim- 
ulations over these durations and combining in the appropriate proportion. The nanoflares 
have a triangular heating profile with a duration of 500 s and amplitude that corresponds to 
a time-averaged energy flux through the footpoints of 5.625 x 10 6 erg cm~ 2 s _1 and 3.0 x 10 7 
erg cm~ 2 s _1 for the weak and strong events, respectively. This time average includes the 
long interval between successive nanoflares. Note that the canonical observed radiative loss 
rate from active regions is 10 7 erg cm~ 2 s _1 (Withbroe & Noyes 1977). In addition to the 
nanoflares, there is a very weak steady background heating that corresponds to a static 
equilibrium with a temperature of 0.45 MK. All strands cool fully (i.e., to 0.45 MK) before 
being reheated by the next nanoflare, but the plasma drains slowly enough that the density 
in the weak nanoflare strands never drops below 15 times the equilibrium density. 

The DEM distribution predicted by the model is given in Figure 2b. The solid curve is 
from the coronal part of the strands, while the dashed curve includes both the corona and 
the transition region footpoints (moss). The agreement with the observed curve in Figure 
2a is remarkably good given the simplicity of the model. It is clear that nanoflare heating 
is consistent with the observations. Note that XRT does not provide significant constraints 
for plasma cooler than about 1 MK (for our filter set), so disagreement with this part of the 
curve is not meaningful. 

Nonequilibrium ionization will reduce the high temperature DEM derived from obser- 
vations relative to the actual values. However, we do not expect the effect to be large for 500 
s nanoflares like those in our simulations (Bradshaw & Cargill 2006; Reale & Orlando 2008). 
The temperature cooling timescale in the strong nanoflare simulation is always greater than 
500 s. 

One can ask whether steady heating is also consistent with the observations. It seems 
likely that a suitable distribution of steady heating could reproduce the DEM curve. How- 
ever, the extreme level of steady heating required to maintain the hottest stands is highly 
implausible. For example, an energy flux of several times 10 9 erg cm -2 s -1 is required to 
produce a static equilibrium at 20 MK. This is more than two orders of magnitude greater 
than the canonical value. If the energy for the heating is provided by the stressing of the 
magnetic field by footpoint motions, then a steady photospheric velocity of order 100 km 
s _1 is required. Such a flow is completely unsupported by observations. Wave heating is 
also implausible. Even with a wave speed as large as 1000 km s^ 1 , the fluctuation velocity 
would need to be of order 1000 km s _1 to provide the necessary energy flux (assuming that 
such waves could dissipate fully in the low corona, which is unlikely). The corresponding 
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displacement amplitude is totally unrealistic for wave periods > 10 s. Thus, we conclude 
that while the cooler plasma in the DEM distribution could be due to either impulsive or 
steady heating, the hotter plasma can only be a result of nanoflares. 

3. Conclusions 

The DEM distribution generated for a quiescent active region core shows a high-temperature 
component with significant hot (Log T > 7.0) emission measure. XRT is sensitive to this 
plasma because it has a greater DEM dynamic range (almost four orders of magnitude) than 
any earlier solar X-ray imager or spectrometer. The hot plasma result is consistent with 
observations from spectrometers in that the synthetic spectrum generated using the XRT 
active region core DEM shows no Fe XIX lines, which would indicate flaring plasma. The 
EBTEL code was used to successfully model the result. For these reasons, we suggest that 
this might be the first convincing evidence of a nanoflare signature. 

This work was inspired by a talk given by Fabio Reale at the September 2008 Hinode- 
2 Meeting in Boulder, CO. We would like to thank Fabio, as well as Paola Testa, John 
Raymond, and Nancy Brickhouse for helpful discussions. Hinode is a Japanese mission 
developed and launched by ISAS/JAXA, with NAOJ as domestic partner and NASA and 
STFC (UK) as international partners. It is operated by these agencies in co-operation 
with ESA and the NSC (Norway). Solar physics research at the University of Memphis is 
supported by a Hinode subcontract from NASA/SAO as well as NSF ATM-0402729. 
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Table I. XRT Data from 2007 May 13 



Filter 


Time 
(UT) 


Long Exposure 

(sec) 


Time 
(UT) 


Short Exposure 
(sec) 


Pred/Obs Flux 
MCMC 


Pred/Obs Flux 
( (xrt_dem_iterative2) ) 


Aljnesh 


18:04:43 


4.10 sec 


18:03:47 


0.18 sec 


1.38±.04 


1.33±.04 


ALpoly 


18:17:06 


5.80 sec 


18:16:38 


0.51 sec 


l 


l 


C-poly 


18:10:00 


8.20 sec 


18:09:01 


0.51 sec 


1.17±.03 


1.12±.03 


Ti-poly 


18:03:27 


8.20 sec 


18:03:17 


0.51 sec 


0.86±.02 


0.82±.02 


ALpoly-Ti 


18:06:33 


16.4 sec 


18:05:23 


1.45 sec 


0.86±.02 


0.81±.02 


C-Ti 


18:08:17 


16.4 sec 


18:07:17 


1.03 sec 


0.80±.02 


0.75±.02 


Bc_thin 


18:12:31 


23.1 sec 


18:11:47 


1.03 sec 


1.47±.04 


1.37±.04 


Be_med 


18:14:06 


46.3 sec 


18:13:13 


2.05 sec 


1.16±.03 


1.09±.03 


Al_med 


18:15:39 


46.3 sec 


18:15:04 


2.05 sec 


1.26±.03 


1.18±.03 


ALthick 


18:02:19 


46.3 sec 


18:01:54 


16.4 sec 


1.99±.07 


1.87±.07 


Be_thick 


18:00:37 


65.5 sec 






1.09±.03 


1.03±.03 



1 Filter not used due to concerns about contamination by atmospheric absorption (near earth eclipse) 
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Fig. 1. — Composite long/short exposure XRT image of AR 10955 made with the TLpoly 
filter. The DEM analysis used contiguous pixels within the outlined area which had >0.03 
Data Numbers s _1 in Be_thick. 
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Fig. 2a. — Optimal median DEM solutions (see text) from xrt_dem_iterative2 (yellow, solid) 
and MCMC (orange, dashed); the minimum \ 2 solutions are similar. The solution sets from 
200 different Monte-Carlo realizations are overplotted as dotted histograms. 
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Fig. 2b. — Theoretical DEM curves generated from a simple composite of two nanoflare 
simulations that have been averaged over time and scaled to match the observations. The 
solid curve includes only the coronal emission where the dashed curve adds the footpoint 
(moss) emission. 
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Fig. 3. — Single-T DEMs with Log T = 6.1 and 6.4 (top) and multi-thermal DEM (bottom). 
The model DEM is shown as a solid line; the grayscale shading indicates the density of 100 
Monte Carlo estimates; and the diamonds indicate the median fit. 
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Fig. 4. — Synthetic spectra generated from CHIANTI using the XRT active region core 
DEM (dot) and the B8 flare DEM (solid) from Reeves et al. ( 2009). The first spectrum is 
dominated by the Ne IX triplet lines at 13.45, 13.55, and 13.70 A as well as an Fe XVII line 
at 13.82 A. The flare spectrum is divided by five for comparison. It shows a strong Fe XIX 
line at 13.52 A as well as other flare lines. 




